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Abstract 

Scaling properties of the cosmic microwave background (CMB) 
radiation are studied using satellite (COBE-DMR maps), balloon- 
borne and ground-based (combined QMASK map) data. Quantita- 
tive consistency is found between the multiscaling properties of the 
COBE-DMR and QMASK CMB maps. Surprisingly, it is found that 
the observed CMB temperature multiscaling resembles quantitatively 
the multiscaling properties of fluid turbulence, that indicates primor- 
dial plasma turbulence as an origin of the CMB temperature space 
anisotropy. 

PACS numbers: 98.70.Vc, 98.80Es, 95.30Qd 
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1 Introduction 



In the last few years the interest in nonlinear (turbulent) processes in pri- 
mordial plasma has been renewed in relation to the origin of magnetic fields 
observed in galaxies (see, for instance |I]-|E] and references therein). It has 
been suggested that primordial magnetic fields might arise during the early 
cosmic phase transitions. The plasma of the early universe has a high con- 
ductivity so that a primordial magnetic field would be imprinted on the 
co-moving plasma and would dissipate very slowly. Such a field could then 
contribute to the seed needed to understand the presently observed galactic 
magnetic fields, which have been measured in both the Milky Way and other 
spiral galaxies, including their halos. At earliest times, the magnetic fields 
are generated by particle physics processes, with length scales typical of par- 
ticle physics. If the inflation hypothesis is correct, long correlation lengths 
can be expected following the inflation. It is shown in [Tj that turbulence 
with its cascade processes is operative, and hence the scale of magnetic fields 
is considerably larger than would be the case if turbulence were ignored. 
The turbulent nature of the magnetic field may have interesting effects on 
various phase transitions in the early universe. Also, the inherent shift of 
energy from small to large scales may be of interest in connection with den- 
sity fluctuations due to the magnetic energy. In particular, it is shown in 
that rotational velocity perturbations, induced by a tangled magnetic 
field can produce significant angular scale anisotropies in cosmic microwave 
background (CMB) radiation through the Doppler effect. The conclusions 
are relevant to arcminute scales [3] and to I > 1000 jlj. These scales will be 
available for observation with the new MAP and Planck space missions. On 
the other hand, the authors of [5] and jH] consider early universe turbulence 
and its imprint on the cosmic background on larger angular scales via a tensor 
mode contribution. The question then is: Can one use the modern data on 
the CMB radiation (see for recent reviews [Zj,[H]) to find out the fingerprints 
of the primordial turbulence? The first attempt to find such fingerprints 
was performed in j^J. However, a CMB map used in Ref. jU] was one of 
the earliest of the COBE-satellite CMB maps and, therefore, the preliminary 
investigation of Ref. P gave a qualitative picture only. In present paper to 
obtain quantitative results we will use the improved four- year COBE-satellite 
CMB maps. Moreover, we have now a new generation of balloon-borne and 
ground-based CMB data obtained with angular resolution of the order of a 
degree (the angle resolution for COBE data, in contrast, was ~ 7°). We 
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choose so-called QMASK map representing the data in a form suitable for 
our purpose. The QMASK map is a combination of QMAP (balloon-borne) 
and Saskatoon (ground-based) data (TUJ- [T3] . Using the three CMB maps 
(two COBE-DMR and one QMASK) we study the structure functions of the 
CMB temperature space increments and moments of CMB dissipation rate 
for frequency ranges: 31.5,53, and 90 GHz (COBE-DMR data), 26-46 GHz 
(QMASK map). 

The main results of our investigation are: 

1) The CMB temperature structure functions and dissipation rate mo- 
ments are found exhibiting multiscaling properties. 

2) The observed multiscaling resembles quantitatively the multiscaling 
properties of fluid turbulence, that indicates its primordial turbulence origin. 

3) Comparison of the results obtained for the COBE-DMR and for the 
QMASK shows quantitative consistency of their CMB maps. 

2 Structure functions of QMASK map 

Let us start from the QMASK CMB map. The QMAP (balloon) and Saska- 
toon (ground-based) data were used to construct this map [TU]- [T3|. The ob- 
servations were made in Ka and Q bands (26-36 GHz and 36-46 GHz respec- 
tively) for different polarization channels. QMAP was designed to measure 
the CMB anisotropy by direct mapping. The Saskatoon data set is different 
in the sense that it does not contain simple sky temperature measurements. 
Instead, it contains different linear combinations of sky temperatures with 
complex set of weighting functions. The QMAP and Saskatoon data were 
combined to produce a CMB map named QMASK. The map was generated 
by subdividing the sky into square pixels of side O ~ 0.3° and consists of 
6495 pixels. The data are represented using three coordinates x, y, z of a 
unit vector R in the direction of the pixel in the map (in equatorial coordi- 
nates). Since some pixels are much noisier than others and there are noise 
correlations between pixels, Wiener filtered maps are more useful than the 
raw ones. Wiener filtering suppresses the noisiest modes in a map and shows 
the signal that is statistically significant. 

The main tool for extracting cosmological information from the CMB 
maps is the angular power spectra. However, such spectra use only a fraction 
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of the information on hand. Taking advantage of the good angular resolution 
of the QMASK map, we can study statistical properties of angular increments 
of CMB temperature for different values of angular separation j!4j . The 
increments are defined as 

AT r = (T(R + r) -T(R)) (1) 

where r is dimensionless vector connecting two pixels of the map separated 
by a distance r, and the structure functions of order p as (|AT r | p ) where (.) 
means a statistical average over the map. 

Figure 1 shows, in log-log scales, the dependence of the moments of differ- 
ent orders, p = 0.2, 0.5, 0.7, 1, 2, 3 of (|AT r | p ) against r for the Wiener-filtered 
QMASK map. The best linear fits are drawn to indicate the scaling 

(\AT r \ p ) ~ r^ p . (2) 

Figure 2 shows the scaling exponents ( p extracted from the QMASK data 
shown in Fig. 1 (circles). We also show in this figure (as dashed curve) the 
scaling exponents ( p corresponding to velocity increments in fluid turbulence. 
For this purpose we use the She-Leveque model |15j . which is in very good 
agreement with the numerous data for velocity increments obtained in inertial 
interval of scales for isotropic fluid turbulence: 

C P = p/9 + 2[l-(2/3) p / 3 ] (3) 

The correspondence seen in Figure 2 is rather remarkable (especially for low 
moments for which the calculations are also the most reliable). 

It appears that for the angular interval between 0.9° and 4° — this being 
the most reliable range for the data ^3] — the temperature increments in 
QMASK map, scale essentially as velocity increments in fluid turbulence. 

It should be noted that under certain circumstances scaling of the mo- 
ments (2) for p = 2 can be related to scaling of corresponding spectrum. 
One of the necessary conditions for such relation is sufficiently fast decrease 
of the spectrum with increasing of wave numbers k > rj" 1 , where r\ is the 
lower edge of the scaling interval of r. For fluid turbulence such fast decay 
of the spectral asymptotic indeed takes place (rj is so-called Kolmogorov or 
viscous scale there) . However, it seems that for the CMB data this is not the 
case and the strong relation between fluid turbulence and the CMB space 
fluctuations takes place in a restricted interval of scales only (and does not 
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determine the spectral asymptotic). 



Background signal can be distinguished from the atmospheric and fore- 
ground contamination signals by their frequency dependence, their frequency 
coherence and their spatial power spectra. Employing these methods, it has 
been shown in a series of papers JO]- JE],5Z|, JE] that all known atmo- 
spheric and foreground sources give a minor contribution to the QMASK 
map. We are thus led to conclude that the background (primordial plasma) 
turbulence is the most plausible reason why the CMB map resembles fluid 
turbulence in its intermittency properties. It should be noted that remark- 
able correspondence between the scaling laws in the near space plasma and 
in isotropic fluid turbulence is a well known puzzle (see, for instance, [T3].[2U] 
and references therein). 



3 Dissipation rate of CMB temperature from 
the four-year COBE-DMR data maps 



The above discovered modulation of the CMB temperature by a fluid turbu- 
lent velocity fluctuations can be explored further. Dissipation rate of kinetic 
energy corresponding to the velocity field u is 



ij 3 

Local space averaging of the dissipation rate is used to construct a space 
dissipation measure |2T] : 



where v r is a subvolume with space-scale r. Scaling law of this measure 
moments, 



(where (...) means an average) is an important characteristic of the dissipa- 
tion rate For isotropic fluid turbulence there exists a simple relation 
between scaling exponents ( p for the velocity increments (see previous sec- 
tion) and the scaling exponents /i p (6) 




(4) 




(6) 



(7) 
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This relationship allows to calculate \i v using the She-Leveque representation 
(3). 

Let us now calculate the dissipation rate of the CMB temperature fluctua- 
tions using the COBE-DMR CMB maps. The COBE satellite was launched 
on 1989 into 900 km altitude sun-synchronous orbit. The Differential Mi- 
crowave Radiometers (DMR) operated for four years of the COBE mission 
and mapped the full sky. The instrument consists of six differential microwave 
radiometers, two nearly independent channels that operate at each of three 
frequencies: 31.5, 53 and 90 GHz. Each differential radiometer measures 
the difference in power received from two directions in the sky separated by 
60°, using a pair of horn antennas. Each antenna has a 7° beam. In this 
investigation we will use a four-year DMR data map of the CMB tempera- 
ture. Most of the Galactic emission has been removed from the map using 
so-called " combination technique" in which a linear combination of the DMR 
maps is used to cancel the Galactic emission [23J. We will also use a DMR 
four- year sky map for frequency 31.5 GHz only. In the last case we will re- 
move the sky regions with violently large temperature fluctuations to remove 
roughly the Galactic contamination (see below and Fig. 4). The first detail 
interpretations of the four- year COBE/DMR data one can find in [2*4]-|28|. 

To organize the DMR data in the temperature maps the sky was divided 
into 6144 equal area pixels. These are formed by constructing a cube with 
each face divided into 32 x 32 = 1024 squares, projected onto a celestial 
sphere in elliptic coordinates. The projection is adjusted to form equal area 
pixels having a solid angle of 47r/6144 sr or 6.7 square degrees. Because the 
7° beamwidth of the sky horn is greater than separation between pixels (2.6° 
average), this binning oversamples the sky. 

Let us introduce a dissipation rate for CMB temperature T in following 

way (cf (4),(5)): 

Xr = — (8) 

v r 

Technically, using the above described pixel map, we will calculate the 
CMB temperature dissipation rate using summation over pixel sets instead 
of integration over subvolumes v r . So that in the multiscaling 

(X P S ) ~ (9) 
the metric scale r is replaced by number of the pixels - s, characterizing size 
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of the summation set (cf (6)). The Xs is a surrogate of the real 3D dissipa- 
tion rate Xr- It is believed that the surrogates can reproduce quantitative 
multiscaling properties of the 3D dissipation rates [22] . 

Figure 3 shows scaling of the CMB temperature dissipation rate moments 
(xf) calculated for the DMR maps (Fig 3a: for combination of frequencies 
31.5, 53 and 90 GHz; Fig. 3b: for frequency 31.5 GHz). Log-log scales are 
chosen in the figure for comparison with scaling equation (9). The straight 
lines (the best fit) are drawn to indicate the scaling. 

Figure 4 shows the scaling exponents \i v extracted from figure 3a (open 
circles) and from figure 3b (crosses), the dashed curve in figure 4 corresponds 
to the intermittency exponents \i v calculated using the She-Leveque model 
[15j for velocity (kinetic energy) fluctuations of fluid turbulence: equations 
(3) and (7). 

Thus, for the satellite data we also observe good agreement with fluid 
turbulence modulation of the CMB temperature. Applicability of the fluid 
turbulence relationship (7) to the CMB data can be itself considered as an 
additional indication of this modulation. Moreover, comparing Figures 2 and 
4 (through the relationship (7)) one can conclude that the QMASK data rep- 
resented in Figure 2 are quantitatively consistent to the COBE-DMR data 
represented in Figure 4. In this respect it should be noted that the QMASK 
and the COBE-DMR data were obtained in about the same range of fre- 
quencies. Very recent suggestions to detect the primordial turbulence using 
the relic gravitational waves 0,1221,1201 might give a complementary exper- 
imental tool for the present observations in the nearest future. 

The author is grateful to the QMAP, Saskatoon, COBE-DMR instru- 
ments teams and to the NASA Goddard Space Flight Center for providing 
the data and support. Discussions on the subject with C.H. Gibson, A. 
Kosowsky and K.R. Sreenivasan as well as the Referee's comments and in- 
formation were useful for this investigation. 
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Figure Captions 



Figure 1. Logarithm of the structure functions (|AT r | p ) against logior for 
the Wiener-filtered QMASK map. The CMB temperature T is measured in 
fiK. The straight lines (the best fit) are drawn to indicate the scaling (2). 

Figure 2. The scaling exponents ( p , corresponding to Eq. (2), for the 
QMASK map (circles). The dashed curve corresponds to the scaling expo- 
nents ( p calculated using the She-Leveque model (3) [T3] for the turbulent 
fluid velocity increments. 

Figure 3. a) The CMB dissipation rate moments (xf) against s (the 
COBE-DMR four-year map for combination of frequencies: 31.5, 53 and 90 
GHz). The dissipation rate Xs is measured in 10 -2 ^K 2 ; b) Analogous data, 
but now for the CMB sky four-year map corresponding to frequency 31.5 
GHz. 

Figure 4. The scaling exponents \i v extracted from figure 3 against p. 
Open circles correspond to the CMB four-year map for combination of fre- 
quencies: 31.5, 53 and 90 GHz, and crosses correspond to the CMB four-year 
sky map for frequency 31.5 GHz. The dashed curve corresponds to the in- 
termittency exponents \i p calculated using the She-Leveque model for the 
kinetic energy dissipation in fluid turbulence (equations (3), (7)) [15j. 



10 



10 7 r 
10 6 - 
10 5 - 
10 4 - 
10 3 - 
10 2 - 

io 1 L 



QMASK (Wiener-filtered) 



p=3 



.01 




p=2 O 



p=1 o- 



o 



.02 



-O 



O 



O 



O 



.03 



o 



-o 



J I I 



.04 .05 .06 .07 



100 r 



QMASK (Wiener-filtered) 



™ " p=0.7 Or 



p=0.5 O 



-o 



o 



p=0.2 O 



o 



o 



o 



-o — o 



J I 



.01 



.02 



.03 



.04 



.05 .06 .07 



■ OQMASK (Wiener-filtered) 
— velocity increments in fluid turbulence 



S 

s 

/ 

/ 

s 

s 



k" 



y 
y 

y 

y 

y 

y 



y\_ 

y 



.5 1 1.5 2 2.5 3 



P 



2 r 



O CM B dissipation (31.5, 53 and 90 GHz) 
V CMB dissipation (31.5 GHz) 
— turbulent energy dissipation 



1.5 - 



1 - 



.5 - 












P 



